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Abstract 

We consider the simplest class of the R-invariant gauge mediation model with the 
gravitino mass in the one to ten keV range. We show that the entropy production 
from the supersymmetry breaking sector makes the gravitino into a warm dark 
matter candidate. We also discuss that the gluino mass can be lighter than the 
wino mass even when the messenger sector satisfies the GUT relations at the GUT 
scale. 



1 Introduction 



Dark matter is an important clue to the theory beyond the standard model (SM). There 
have been proposed a lot of models to explain dark matter. The light gravitino is a 
very interesting candidate for dark matter among them, since the gravitino itself is a 
unique and inevitable prediction of supergravity (SUGRA). If the gravitinos were in the 
thermal equilibrium in the early universe, the gravitino mass m^i2 is required to be ^3/2 — 
100 eV from the observed dark matter density, ^^dm — 0.1. This prediction is very much 
interesting, since we can test the gravitino dark matter hypothesis at LHC The gravitino 
mass, m3/2 — 100 eV, is, however, too small to be the cold dark matter and is disfavored 
for the successful galaxy formation [1]. 

The above argument is based on an unjustified assumption on the thermal history of 
the early universe. In fact, if we had late time entropy production after the decoupling 
time of the gravitino, the mass of the gravitino dark matter may be raised up to a few 
keV. Moreover, the gravitino dark matter with a mass in the one to ten keV range serves 
as the warm dark matter which has recently been invoked as possible solutions to the 
seeming discrepancies between the observation and the simulated results of the galaxy 
formation based on the cold dark matter scenario [2] El 

In this paper, we discuss the late-time entropy production from the SUSY breaking 
sector. As we will show, large entropy can be produced from the SUSY breaking sector 
when the sector has meta-stable particles whose lifetimes are long enough to dominate 
the energy density of the universe before they decay. As a result, the gravitino with a 
mass in the one to ten keV range can be a good candidate for the warm dark matter with 
the help of the entropy production. 

The gravitino mass in a range of the one to ten keV also has an interesting impli- 
cations on the phenomenological aspects of the supersymmetric standard model (SSM). 
For the gravitino mass in the keV range, we are led to consider the models with gauge 
mediation [5], E] where the SUSY breaking effects are mediated to the SSM sector via the 
gauge interactions. 

^ The detailed analyses in Refs. [31 14] have placed lower bounds on the warm dark matter mass around 
a few keV range. Thus, it is safer to assume that the dark matter is in the ten keV range. 
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In this paper, as a particular example of the models with gauge mediation, we consider 
a class of the direct mediation models developed in Ref. |7], where the SUSY breaking 
vacuum is stable. The important feature of this class of models is that the models possess 
an R-symmetry. It should be noted that, independent of the SUSY-breaking mediation 
scheme, the (discrete) R-symmetry is considered to be a crucial symmetry for any low- 
energy SUSY extension of the Standard Model. This can be seen from the fact that SUSY 
should be broken at very high energy scale to obtain the nearly vanishing cosmological 
constant if the R-symmetry is largely broken by the constant term of the superpotential 
in supergravity. Therefore, it is quite tempting to consider a mediation mechanism which 
possesses an R-symmetry. 

The notable feature of this class of gauge mediation models is a peculiar spectrum of 
the superparticles. Especially, the gaugino masses do not satisfy the so-called the Grand 
Unified Theory (GUT) relations even if the masses and the couplings of the messenger 
fields satisfying the GUT relations at the GUT scale f8]. For example, the light gluino 
of mass 300 GeV— 1 TeV is achieved with the heavier wino of mass 500 GeV— 2 TeV even 
for the boundary condition satisfying the GUT relations at the GUT scale. Such a light 
gluino will be easily produced at the LHC, and hence, almost all the parameter space is 
expected to be probed by the LHC experiment. 

The organization of the paper is as follows. In section^ we discuss the mass spectrum 
of the SSM particles in the simplest class of the R-invariant gauge mediation model for 
the gravitino with a mass in the one to ten keV range. There, we show that the typical 
gaugino mass spectrum is distinctive from the so-called minimal gauge mediation model. 
In section[3l we discuss an entropy production mechanism which makes the gravitino dark 
matter scenario with a mass in the one to ten keV range consistent with the observed 
dark matter density. The final section is devoted to our conclusions. 

2 An R-invariant gauge mediation model 

Let us discuss the minimal R-invariant gauge mediation model developed in Ref. [7J. We 
introduce two pairs of massive messengers, \E'i, with i = 1,2. Here, and \E'i transform 
as 5 and 5* in terms of the minimal SU{5) GUT representations, respectively. We further 
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introduce a SUSY-breaking gauge singlet field S which has non-vanishing expectation 
values of the F and A terms, 

{S{x,9)) = {S)+F9'. (1) 

(We abuse the notation for chiral fields and its lowest components.) We assume, through- 
out this paper, that the F term is the dominant component of the SUSY breaking and 
hence the gravitino mass is given by 

"''^ " tSt, 

Here, Mp ~ 2.4 x 10^^ GeV is the reduced Planck mass. 

Let us assume that the superpotential in the messenger sector is given by, 

kS m \ / \l'i 



where k and m denote the coupling constant and the mass parameter, respectively. We 
see that the above superpotential is invariant under an R-symmetry with the charge 
assignment, S{2), \E'i(0), \^i(0), ^'2(2), ^^2(2). Notice that the vacuum expectation 
value of the scalar component of 5* breaks the R-symmetry spontaneously. 

So far, we have treated messenger fields in an SU (5) gvt symmetric way. Below the 
GUT scale, however, the 5't/(5)GUT messenger multiplets split into \E'j — )■ and 

-)■ which transform as (3_i/3, 21/2) and (3i/3,2_i/2) under the SM gauge 

groups, SU (3)c X SU (2)l x f/ (l)y , respectively. In the followings, we name the messengers 
with the superscripts d and i "down-type" and "lepton-type" , respectively. In accordance 
with the above splitting, the coupling constants and the mass parameters in Eq. ([3]) may 
take different values for each type of messengers, i.e. k^'^'^^ and m^'^'^\ Especially, the 
renormalization group (RG) evolution makes them different at the lower energy scale, 
even if we impose k^^^ = k^^^ and m^^^ = m^^^ at the GUT scale. 

In this model, we can take k^^\ m^^\ (S) and F as real positive by the phase rotation 
of the fields without loss of generality. Therefore, we can avoid CP violation in the present 
model. In the following of this paper, we take these parameters as real positive. 

Let us discuss the condition for the messenger scalar not to be tachyonic. The mass 
parameters are required to satisfy k^'^^ F / m^'^^'^ < 1 and k^^^ F / m^^^'^ < 1 for the messenger 
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fields not to be tachyonic. By the RG effect, tfie condition for tfie lepton-type messenger 
gives severer constraints tfian tfie one for tlie down-type messenger. Tliis fact can be seen 
as follows. The ratio of k^'^^ F / m^'^^^ and k^^^F/m^^^"^ evolves according the RG equations, 

d fk^'^^F /k^^^F\ 



~ l^^_3^_i^ (A) 
3 Att Att 3 An' ^ ' 

where we have used the RG equations of k and m in terms of the anomalous dimensions 
of *S^)(^l^Vnd 5, 7*(x) and 75, 

i 

^ A;W = (2^^(„+^5)A;W, (5) 



9 log n 



d 



-m 



- (7^(x)+7^(x))mW, (x = c?,£). (6) 



9 log IX 

If we require = fc*^^^ and mS'^^ = m^^^ at the GUT scale, we can get 



Mgut / 16 as 1 ai 

a log u : 6- - — 

M^ed V 3 47r 47r 347r 



(7) 



at the mediation scale. Due to the strong SlJ{3)c effect, k^'^^ F j m^^^'^ becomes smaller than 

k^^^ F I vrS^'^'^ at the mediation scale. Therefore, the condition k^^"^ F j vn!^^^'^ < 1 guarantees 
that all the messenger scalars have positive squared masses. The fact k'-'^^ F / m^'^^^ < 
/jW^y/^(^)2 -j^Qg Q^j^ interesting consequence for the gaugino masses as will see later. 

Mass spectrum of the present model 

The gaugino masses are given by 

M. - ^(fAS?.4AS^). («) 

= (9) 
Af3 = gA<'A, (10) 

and the squared masses of sfermion / are given by 



^~ Ka-k) 5 
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where aa [a = 1,2,3) are gauge coupling fine structure constants of f/(l)y, SU{2)l, 
SU{3)c, and Ca {a = 1,2,3) are quadratic casimir invariantj^ of the sfermion / under 



and A^^^' 

m^^\ /^(x) ^5*^ g^Yid k(x)p whose exphcit forms can be read from Ref. [HI E]. 



the group U{1)y, SU{2)l, SU{3)c. Here, A^)^ and A^"^^" (x = dj) are functions of 



Gaugino-sfermion mass ratio 

In the case of the so-called minimal gauge mediation (mGM) [6], with a messenger super- 
potential of the form (m*^^^ + k'^^^S)'^'^, both the A^^\ and Aq^"* are of order 



A 



(x) I 

1/2 I mGM 



Ao'^'*|mGM ~ 



kix)F 



m 



(x) 



(12) 



where we have neglected (S) (inclusion of it is straightforward). However, there is a 
significant difference in the gaugino masses in the present model. The rough behavior 
of the soft masses are as follows. For A;*-^^ (S) ^ m*--^-*, the soft masses aS)^ and aJ^)' are 
approximately given by 



A 



(x) 

1/2 



O(O.l) X 



(fcW (5))(A;WF) 



A, 



(X)2 



0(1) X 



( 




(13) 



(14) 



and they decrease when k^^'^S becomes much larger than m*--^-*. Notice that there are no 
terms of order k^^^F/m^^^ in the gaugino masses, which would be present in the minimal 
gauge mediation. Because the mass parameters are required to satisfy k^^^F < m^^^'^ for 
the messenger fields not to be tachyonic, the above approximated expressions show that 
the gaugino masses are suppressed compared with the sfermion masses. 

Wino-gluino mass ratio 

In the case of the minimal gauge mediation, the ratio of gaugino mass contributions from 
the down-type and the lepton-type messengers, A['^2lmGM and A^^^gUcM, is given by 



^^1/2 



m 



m 



(15) 



mGM 



We use the GUT normalization for the U{1)y gauge group. In particular, the quadratic casimir is 

5' 



given by Ci = ^Y^ in terms of the hypercharge Y. 
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where we have neglected the higher order terms in k^^^ F / m^^'^'^ . One can easily check 
that this ratio is invariant under the RG flo\\c|. Thus, if we impose the GUT relations, 
f^{d) _ f^{e) g^^^ ^{d) _ ^(t) ^ GUT scale, we have A^y2 — ^1/2 messenger scale. 

This indicates that the gaugino masses obey the famous GUT relation Mbmo : M^ino : 
MgMno ^ ai : a2 : as - 1 : 2 : 6. 

In the present model, on the other hand, the ratio A^y2/^i/2 given by 



(16) 



By using similar argument to derive Eq. ([7]), the ratio Eq. ( !T6|) at the scale of the gauge 
mediation, M^aed = 0{m^'^'^^), is roughly given by 

Ad) 



exp 



exp 



/ rflog/i 2(7^(d) - 7^(.)) - 6(7^{d) - 7^w) 
-AfcuT ^ ^ / 32 0(3 g 2 



;i7) 



if we impose the GUT relations at the GUT scale, Mqut, i-e. m*^'^) = m*^^^ and /e*^'^'' = k'^^\ 
The strong 5'f/(3)c interaction makes A^^g smaller than A[^y*2 the mediation scale. 
Therefore, from Eqs. fl8|)-f lT0|) . we see that the gluino mass is rather suppressed [8] in 
the present model compared with the one in the minimal gauge mediation. Furthermore, 
when fc^F ~ m^^^^, the ratio AS/yAS/^2 

is more suppressed than Eq. ( ITTl) . This is because 
the higher order contributions of fc'^^-'F/m^^)^ in Eq. f|T3|) make A[^2 enhanced, while A[^2 
is not so enhanced since the RG equation indicates k^'^^ F / m^'^^'^ < k^^^ F / m^^^'^ as we have 
discussed above. 

Numerical results 

We show numerical results of the soft SUSY breaking masses of gauginos and sfermions in 
the minimal R-invariant gauge mediation model. We impose m^'^^ = m^^^ and k^^^ = k^^^ 
at the GUT scale. The larger soft masses are obtained for the larger coupling constants. 



^ As pointed out in Ref. [TU], this RG invariance is a general feature of models where gaugino masses 
are generated at the leading order in the SUSY breaking F-term. In the effective field theory language, 
the leading term in F is generated [TT] by a holomorphic gauge kinetic function J (P9h{S)W°'Wa and 
this holomorphic term is independent of the wave function renormalization of the messenger fields. 
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The soft SUSY breaking masses of gauginos and sfermions in the minimal R- 
model for 1713/2 = IkeV (left) and ms/2 = lOkeV (right). We set k^^^ F / m^^^'^ = 
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Figure 2: The soft SUSY breaking masses of gaugino and sfermion in the minimal R- 
invariant model for 7713/2 = IkeV (left) and 7/73/2 = lOkeV (right). We take (S) as a value 
which maximizes the gluino mass. 
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k^"^^ and fc'-^^ at the mediation scale. The couphng constants are, however, not able to be 
arbitrary large, because the Yukawa-type interactions are not asymptotically free. Too 
large coupling constants result in the Landau-pole problem below the GUT scale. In 
order to avoid the Landau-pole problem, we put k^^^ = fc*-^-* = An at the GUT scale as 
the upper bound on With this boundary condition, we obtain 

= 0.99, k^^^ = 0.74 and m^'^^/m'^^^ = 0.71 at the mediation scale in the minimal R- 
invariant gauge mediation model for m3/2 = — (9(10) kev|^ The explicit form of the 
RG equations are given in Ref. [9]. In this analysis, we have neglected the contribution 
of the SUSY breaking sector to the RG equations, which may make the values of k^^^ a 
little smaller. 

In Fig.[T], we show the SSM mass spectrum as a function of the parameter {S) /mS^'^ for 
^3/2 = 1 keV and ^3/2 = 10 keV. We fixed k^^^F/m^^^'^ = 0.9 in both figures. The figures 
show that the sfermion masses are much heavier than the gauginos. The sfermion-gaugino 
mass ratio is larger for the heavier gravitino mass if we fix the order of the gaugino masses^ 
The figure also shows that the gaugino masses are maximized when k^^^S) ~ m^^\ In 
Fig. [21 we also show the mass spectrum of gaugino and sfermion as a function of the 
parameter k^^^ F / m^^'^'^ . In the figure, we took {S) as a value which maximizes the gluino 
mass. 

The interesting observation here is that the gluino mass is predicted to be rather 
light and can be lighter than the wino even though we have assumed the GUT boundary 
conditions. Thus, the gaugino mass spectrum is distinguishable from the one in the 
minimal gauge mediation, Mbino : ^wino : Mgi^i^o — 1:2:6. 

In Fig.OandlU we also show the SSM mass spectrum for the model with two additional 
pairs of massive messengers \Ef^, "^'^ {i = 1,2) which have the masses and the couplings 
to 5* similar to those of \E'j and \E'j in Eq. ([3]). In the next section, we consider the SUSY 
breaking model where the R-symmetry is spontaneously broken in a perturbative way with 
the help of the U{1) gauge interaction. In that model, the messenger sector is required 

^ The perturbative analysis is no more viable for k — 0{A7t). Our result, however, does not strongly 
depend on the values of the coupling constants at the GUT scale as long as they are large. 

^ As we see shortly, the mediation scale is required to be close to y/F to obtain heavy enough gaugino 
masses, which is determined for a given gravitino mass. 

^ In Refs.[T31 [H], the SSM spectrum in the R- invariant gauge mediation with the very light gravitino 
with a mass below 16 eV has been considered. There, the sfermion-gaugino mass ratio is much smaller. 
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Figure 3: The soft SUSY breaking masses of gauginos and sfermions in the R-invariant 
model with double messenger for m^i2 = IkeV (left) and 171^/2 = lOkeV (right). We set 
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Figure 4: The soft SUSY breaking masses of gaugino and sfermion in the R-invariant 
model with double messenger for m3/2 = IkeV (left) and = lOkeV (right). We take 
(S) as a value which maximizes the gluino mass. 
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to be doubled (see also Eq. (!70l) ). The sfermion-gaugino mass ratio is smaller than the 
minimal model, because the gaugino mass is proportional to the messenger flavor number 
but the sfermion mass is proportional to square root of the messenger flavor number. 

The flgure shows that the gluino mass is more close to the bino mass in the doubled 
messenger model. For such a peculiar spectrum, the collider phenomenology can be 
signiflcantly different from the usual models with the gaugino masses which satisfles GUT 
relations, and hence, may require different search strategies than the usual SUSY scenarios 
(see for example Ref. [12] )• The detailed collider study for the above gaugino spectrum 
will be given elsewhere [13] . 

3 Entropy production from SUSY breaking sector 

In the previous section, we have discussed the mass spectrum in the minimal R-invariant 
gauge mediation model, and the peculiar mass spectrum are predicted for the gravitino 
mass in the one to ten keV range. As we have mentioned in the introduction, however, 
the relic density of the gravitino with a mass in this range is too high to be consistent 
with the observed dark matter density if the gravitinos were in the thermal equilibrium 
in the early universe. The thermally produced gravitino density is roughly given by, 

\9*iTD)) VlOOeVy ' ^ ^ 

where g*{T£,) ~ 100 denotes the effective massless degree of freedom in the thermal bath 
at the decoupling temperature. To, of the gravitino from the thermal bath [16], 

1/2 



To ~ max 



M_.aOev(^)'(i^)X^)- 



(19) 



Notice that the effective interactions between the gravitino and the SM fermions after 
integrating the SUSY particles out are so suppressed that they cannot keep the gravitino 
in the thermal bath after the gauginos decouple. In the above expressions, we have 
neglected the contributions from the Winos and Binos which could give comparable or 
even a larger contributions. The following discussion is not affected as long as the order 
of the magnitude of To is not changedlll 



One may obtain more precise expressions of Tjj by using more recent analysis on the gravitino 
production cross section given in Ref. [17) . 
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From the above discussion, for the gravitino with a mass in the one to ten keV range 
to be a consistent dark matter candidate, the above rehc density should be diluted by 



after the decoupling of the gravitinojj It should be noted that the late time entropy 
production dilutes also the primordial baryon asymmetry by the same factor, but it may 
not cause any serious problem [TS] in the thermal leptogenesis for A ~ 10 — 100. 

In the previous section, we have also made a tacit but a crucial assumption; the 
spontaneous R-symmetry breaking, (S) ^ 0. It is not trivially realized in many SUSY 
breaking models. 

In this section, we propose a very ambitious solution to both the above problems, 
the dilution of the thermal gravitino and the R-symmetry breaking, at the same time, by 
considering the entropy production from the SUSY breaking sector, where the spontaneous 
R-symmetry breaking is realized. 

3.1 Extended vector-like SUSY breaking sector 

As an example of the SUSY breaking model where the spontaneous R-symmetry breaking 
is realized, we consider the extended model of the vector-like dynamical SUSY breaking 
model based on SU{2) gauge theory in Ref. [2T|, |22]. In the extended model, one of the 
global U{1) symmetry is upgraded to a gauge symmetry [23j which is spontaneously broken 
at the SUSY breaking vacuum. As discussed in Ref. [23], the spontaneous R-symmetry 
breaking is achieved in a perturbative way, which is not the case in the original SUSY 
breaking model. 

The notable property of the extended model is that it possesses an accidental discrete 
symmetry even after the spontaneous U{1) symmetry breaking |24j . Thus, the lightest 
particle which is charged under the unbroken discrete symmetry has a long lifetime. As 
we will show shortly, the energy density of such a long lived particle can dominate the 
universe and cause the entropy production when it decays, which dilutes the thermally 
produced gravitino. 

^ The gauge mediation mechanism often involves natural mechanisms of late time entropy production 
from, for example, the messenger sector [18^ or the intermediate SUSY breaking sector [19j. 




(20) 
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Table 1: The symmetries of the model. SU (2) x f/(l) are gauge symmetries and Z4 x f/(l)H 
are global symmetries. Notice that both the global symmetries are anomaly free. 





SU{2) 


f/(l) 


^4 




'S'12 


1 


-1 




2 


5'34 


1 


1 




2 


'S'la, 14,23, 24 


1 







2 


Ql,2 


2 


1/2 


gi7r/2 





^3,4 


2 


-1/2 


gj7r/2 






The vector-like SUSY breaking model consists of four fundamental representations 
Qk{k = 1, ■ ■ ■ , 4) and six singlets Sij = —Sji = 1, ■ ■ ■ , 4) which interact with each 
other in the superpotential, 

W = Y.X'^^S,,QkQi , (21) 

where A's denote the coupling constants. 

In the extended model, we gauge one of the U{1) subgroup of the maximal subgroup 
of SU{4:) of the model by assigning the gauge charges, Qi^2(l/2), Q3^4,{— 1/2), S'i2(— 1), 
5'34(+l), and 5*13,14,23,24(0). With this charge assignment, the above superpotential reduces 
to 

w = Aw^i2gig2 + x'^-^SuQsQ, + E ^1is^JQkQl , (22) 

where A'(^j = for ij = 12, 34 or kl = 12, 34. The global symmetries of the model are 
f/(l)/j X Z4. The charge assignment of the R-symmetry is S{2) and Q{0). The fields are 
transformed to iQ and —5* under the Z4 symmetry]^ We list the symmetries of the model 
in Tabled] 

Below the dynamical scale A, the model is well described by using the composite fields, 
Mij ~ QiQj/A, whose superpotential terms are approximated by, 

Wes = A(+)A^i2Mi2 + A(-)A^34M34 + E >^'ij^SijMki + A'(Pf (M,,) - A^) , (23) 

^Classically, the Z4 symmetry can be realized as a continuous U{1) symmetry with the charge as- 
signment S{2) and Q{—1). The anomaly against the SU{2) gauge symmetry breaks the U{1) symmetry 
down to the discrete Z4 subgroup. 
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where X denotes the Lagrange multipher which expresses the quantum deformed moduh 
constraint, Pf(M) = A^. Here, the ambiguity of the normahzations of the meson fields 
are imphcitly absorbed by A's. Thus, strictly speaking, the parameters A's appearing in 
Eqs. ( |23ll and (12^ are different from the ones in Eq. (1221) . By using appropriate linear 
combinations of 5"s and M's we may rewrite the above effective field theory by. 



W = A(+)AS+M_ + A(-)A5_M+ + E K^SaMa + X M+M_ + ^ ^MaM,, - A^ .(24) 

Here, the newly introduced matrix yah is generically given by, 

= iJJ^V)ab. U e SU{A) . (25) 

By assuming that A's are perturbative, and A^ are smaller than A"s, we may parametrize 
the deformed moduli space, byjl^ 



A2- E —MaMb, M_=e-^/^^ 

a,b=1...4 2 



\ 



A^- E ^AifoMj. (26) 



Notice that the U{1) gauge symmetry is spontaneously broken on the deformed moduli 
space. Then, the above effective theory can be reduced to 



W = (A(+)A5+e-'^/^^ + A(-)A5_e<^/^^) 



A2 _ ^ Mm,M, + E K^SaMa .(27) 

a,b=1...4 ^ a=l-4 



In the foUowings, we simplify the model by taking A*^^^ = A and yat = Sab, although we 
can easily generalize the results for more generic coupling constants. 

By expanding the above superpotential around = and Ma = 0, we obtain. 



E mA + E K^SaMa 
a=l---4 / a=l---4 



v^AA^S - AAT0 + -^S U - E Ml] + E K^SaMa 



{2i 



where we have defined S and T by, 

1 



1 



S = —{S+ + S^), T=-^iS^-S^ 



(29) 



-'^"The "radial" component of the M± becomes a "mass partner" of X, and hence, we can mtegrate the 
radial component out. 
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Table 2: The symmetries of the model in terms of the low energy fields. By the conden- 
sation, (M_(.M_) = A^, the symmetries U{1) x Z4 break down to a global Z2 symmetry 
under which S„ and M„ are odd. 





f/(l) 


z. 


U{1)r 




-1 




2 




1 




2 









2 


M_ 


-1 


giTT 





M+ 


1 


















The newly defined S corresponds to the pseudo-flat direction which F-term breaks SUSY. 

The above effective theory possesses an discrete symmetry, Z2, after the SUSY and the 
f/(l) gauge symmetry breaking, which is a diagonal subgroup of the discrete subgroup of 
the U{1) gauge symmetry and the global Z4 symmetry. Under the Z2 symmetry. Ma and 
Sa are odd while the other T, and are even (see Table[2|)J^^I Therefore, the lightest 
particle which is odd under Z2 symmetry is stable. In the appendix, we give the mass 
spectrum of the SUSY breaking sector which shows that the lightest Z2 odd particle is an 
appropriate linear combination of the scalar components of Ma and 5*^. In the foUowings, 
we name the lightest Z2 odd particle, x^. 

3.2 Decay of the Z2 charged particle 

So far, we have assumed that the Z2 accidental symmetry is an exact symmetry. Generi- 
cally, such an accidental symmetry could be broken by higher dimensional operators which 
are suppressed by the reduced Planck scale. For example, the lowest dimensional operator 
which breaks the Z2 symmetry is 

c cA 
W ~ —QQH^H, ~ MH^H, , (30) 

IVlp IVIp 



In addition to the Z2 symmetry, the model apparently possesses another discrete symmetry under 
which (f) and T are odd. This symmetry stems from the charge conjugation symmetry of the gauged U{1) 
symmetry. The charge conjugation symmetry is, however, expected to be broken in generic models with 
Xi+) ^ xi-)^ In the following s, we do not consider the charge conjugation symmetry as a good symmetry. 
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Figure 5: Left) The yield of the hghtest Zi odd particle, xl, for a given SUSY breaking 
scale. Right) The domination temperature of the lightest charged particle after freeze 
out. In the figure, we have used g^, ~ 250 at the time of the freeze-out of xl- 

where c is a coupling constant and Hu and Hd denote the Higgs doublets in the SSM. 
Here, we have assumed that the i?-charge of HuHd is 2, assuming that the R-symmetry 
(or its discrete subgroup) is better symmetry than the accidental Z2 symmetry. 

Through this operator, the lightest Z2 odd particle decays into Higgs (Higgsino) with 
the decay rate 

1 c^A^ 

Ta-r ~ TT'^nrixT ■ (31) 

Here, we have neglected the masses of the Higgs and Higgsinos in the final state. As a 
result, the decay temperature is roughly given byl^ 

...^(^)-'*ys:5;^icev..0'7^)(^)"'^ . (.) 

Therefore, the decay temperature of the lightest Z2 odd particle is expected to be very 
low even if the Z2 symmetry is explicitly broken by the higher dimensional operators. 



3.3 The entropy production and dilution of the gravitino 

Since the decay temperature of the lightest Z2 odd particle is very low, the energy density 
of the lightest Z2 odd particle is expected to dominate the universe, where the domination 
temperature is estimated by, 

4 



^The detailed analysis is given in the appendix. 
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Figure 6: The dilution factor for a given SUSY breaking scale. In the figure, we have 
fixed Tdecay — 100 Me V by choosing an appropriate coefficient c in Eq. (15^ . 

where is the yield (number density divided by entropy density) of xl after its freeze- 
out. The yield is roughly given by, 



Y^^ ~ min 



0.278 0.76 



(34) 



9* gV^MpTf (awrei) 

where g^, is the effective massless degree of freedom at the freeze-out temperature Tj, and 
the (cTfrci) is the thermal averaged annihilation cross section, 

1 A"^ 

(C^^^rel) ~ T, — ■ (35) 

As we show in the appendix, the dominant mode of the annihilation process of is the 
one into a pair of the gravitinosl^ 

In Fig.|5l we show the yield and the domination temperature. The figure shows that 
the domination temperature is typically below the decoupling temperature of the gravitino 
given in Eq. f ll9p . Therefore, after the decay of xl, the yield of the gravitino is diluted by 
the factor A~^, 

As after decay 4 p after decay 4 TTlx j Tlx r before decay 4 
= —r — ~ r — ~ , ~ rrirrY^r 

■5 1 before decay "J ^ decay ^ \ before decay "J ^ decay ^ \ before decay "J decay 

where s is the entropy density, p is the radiation energy density, and rtx^ is the number 



"'^^The gravitinos produced at the annihilation of interact with the thermal bath and are thermalized 
immediately since T/ ^ To- 



ll 



density of xl- Thus we obtain (using Eq. (|33ll ) 

A ~ . (36) 

-^decay 

In Fig.El we show the resultant dilution factor for a given SUSY breaking scale. In the 
figure, we assumed Tdccay — 100 MeV which roughly corresponds to the lower bound on 
the decay temperature not to affect the Big-Bang Nucleosynthesis, 3 The figure shows 
that the dilution factor of the order of 10 — 100 is achieved for A ~ 0(10^^), which is 
required to achieve the consistent gravitino dark matter scenario with a mass in the one 
to ten keV range. 

3.4 Gravitino from the process of entropy production 

The thermal relic gravitinos are diluted as we have discussed above. In this subsection, we 
discuss other sources of gravitinos. The decay of xl produces higgsinos with a branching 
ratio of order 1, and they eventually decay into gravitinos after some cascade decay. 
Besides, the particles from the xl decay have very large energies of order rrixj^, and they 
interact with the thermal bath and may produce more SUSY particles. 

Suppose that A^^;^ gravitinos are produced in average from a decay of single xl. Then 
the yield of the gravitinos coming from the decay is, 

-i^decay t\t ^xj^ -^xl Plafter decay 3 ^decay /o'?\ 

^3/2 ~ -\ -— -\ jN^L — — • (37) 

* I after decay '"'Zi "'latter decay ^ '"'Xi 

Thus, the contribution of these gravitinos to Qh'^ is, 

^ VIO keVAlOO MeVy V rn^L J 
The decay remnants of xl whose initial energies are of the order of -E ~ rrixj^ lose most 
of their energies by interacting with the particles in the thermal bath. Then, the SUSY 
particle productions take place when the energies of the remnants decrease down to the 
threshold energy, E ^ ^gaugino /^decay As a result, the average number of the gravitinos 
from the decay of is at most of the order of unity [18] . 

Here, instead of analyzing the detail of the above process, we give a very rough upper 
bound on N^j^ by neglecting the energy loss of the remnants, which is good enough for 



The decay temperature, Tdocay — 100 MeV, can be achieved by choosing appropriate coefficient c. 
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our discussion. That is, the maximum number of the SUSY particles from a decay of 
xl is achieved if the initial energy E ^ rrix,^ is distributed to iV ~ "^xi/l^gaugino/^decay) 
particles with energy E ~ ^gaugino/^decay, and each particle with this energy produces a 
gaugino. Thus, the absolute upper bound on the number N^^ is given by, 

N..<N ^ lo. AOOGeV V/ T.eeay W m V 2Br(a:,^i^„g,), (39) 
V ^gaugino/ VlOO MeVAlOe GeVy ^ u d). \ i 

where the second term comes from the direct higgsino production by the decay oi xl- 
Therefore, we find that the gravitino coming from the xl decay does not become a dom- 
inant component of the dark matterl^ 

3.5 Spontaneous R-symmetry breaking 

Finally, we discuss the R-symmetry breaking in the extended vector-like SUSY breaking 
model. As we see from the superpotential in Eq. (1281) . the SUSY breaking field S has the 
flat potential at the tree level, and hence, S corresponds to the pseudo-flat direction. The 
potential of 5* is, however, deformed by the radiative correction, and especially, the origin 
of S can be destabilized |23j . 

At the one-loop level, the effective potential of the pseudo-fiat direction is given by 
the so-called Coleman- Weinberg potential, 

Vc^{S) = -— tT{-fM\S) log , (40) 

where fiR denotes the renormalization scale, and {—)^ = 1 for bosons and (— )^ = —1 for 
fermions. The mass spectrum in the extended model is given in the appendix. 

In Fig. [71 we show the parameter space where the spontaneous R-symmetry breaking 
is achieved for given values of A. In the figure, gx denotes the gauge coupling constant of 
the U{1) gauge interaction. The figure shows that spontaneous R-symmetry breaking is 
realized for gx ^ A, k^^\ which corresponds to the large contribution from the U{1) gauge 
interaction. Notice that for (7x~0.5 at the mediation scale, the U{1) gauge interaction 
has the Landau pole below the GUT scale. Thus, if we require that the extended model 



"'^^ Although the energy density of the gravitino component from the decay of xl is subdominant, such 
gravitinos may have much larger velocity than the one of the thermally produced gravitino. Thus, the 
gravitino component from the decay of xl could have some impacts on the structure formation. 
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is perturbative up to the GUT scale, the gauge couphng constant should satisfy gx ~ 0.5 
at the mediation scale. 

For the larger field value of S* 3> Mmed, the radiatively generated potential is further 
approximated by [25j, 



V^cw(^) 



IF 



2 



Zs 

Z{S) = exp 



In 5 

275 dhifj, 

In A/mcd 



(41) 



Here, 75- is the anomalous dimension of 5* which is, at the one-loop level, given by, 

275 = + + 0{X') , 42 

where (9(A^) contribution comes from mesons M in the region XS ^ A or quarks Q in the 
region AS* ^ A, and this contribution is always positive in the region where perturbative 
calculation is valid (i.e. AS* is not near A). In order for the potential not to show the 
runaway behavior, we need to have 75 > at the larger value of S. Because (9(A^) 
contribution is positive, we can conservatively neglect this contribution to constrain the 
viable parameter space. So we neglect this contribution in 75. In Fig.[71 we have shown 
the region where 75 > at least around the mediation scale. The figure shows that 
spontaneous R-breaking is realized for g v ^ A, k'-^\ while the runaway behavior is avoided 
for not too large gx compared with 

Put the above discussions together, we find that spontaneous R-symmetry breaking 
and the right amount of the entropy production are achieved at the same time for A = 
O(10~^) in the extended vector-like SUSY breaking model. 

4 Conclusion 

In this paper, we discussed the simplest class of the R-invariant gauge mediation model 
for the gravitino with a mass in the one to ten keV range. The gravitino dark matter 



In the figure, we have fixed V k^^'iF jra^^^ — 0.8. For the heavier messenger masses, the hatched 
regions in the figure are shifted to right while the light-shaded region is not shifted. This is because 
the contribution to the Coleman- Weinberg potential in Eq. (j40p from the messenger is suppressed by the 
heavier masses. Thus, for the heavier messengers, the allowed parameter space is larger. 
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Figure 7: Spontaneous R-symmetry breaking is realized by the radiatively generated 
potential of S in the hatched region for given values of A. The U{1) gauge coupling gx 
and the messenger- SUSY breaking field coupling k^^^ are defined at the mediation scale. 
In the figure, we have fixed A' = 2A, k^"^^ ~ 1.8k^^\ m^'^) = 2.1m(^), and = 
We have fixed the value of F. In the dark shaded region, the U{1) gauge coupling constant 
at the mediation scale is too strong and has a Landau pole problem below the GUT scale, 
i.e. gx ~ 0.5. In the light-shaded region, the radiatively generated potential does not curl 
up for the large value of S. 

scenario with the mass in this range is drawing attention as an interesting interpretation 
of seeming discrepancies between the observation and the simulation of the structure 
formation based on the cold dark matter model. 

For a consistent gravitino dark matter scenario with a mass in the one to ten keV range, 
the relic density of the gravitino is needed to be diluted by a factor of A = 10 — 100. In 
this paper, we discussed entropy production from the vector-like SUSY breaking model, 
which is extended so that spontaneous R-symmetry breaking is achieved. Spontaneous R- 
symmetry breaking is necessary ingredient for the successful R-invariant gauge mediation 
mechanism. As a result, we find that R-symmetry breaking and right amount of entropy 
production are achieved at the same time for a certain parameter space. 

The interesting prediction of the R-invariant gauge mediation model is the peculiar 
gaugino mass spectrum with much heavier sfermions. Especially, the gluino can be lighter 
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than the wino even if the messenger masses and couphng constants satisfy the GUT rela- 
tion at the GUT scale. The light gluino (mgiuino — 300 GeV — 1 TeV) is quite advantageous 
to be produce at the LHC. Therefore, it is expected that this model can be probed by the 
LHC. 

Another notable point is that the solution to the yU-problem proposed in Ref. [26] 
works for the gravitino mass of order 0{1) — (9(10) keV. It is very interesting that the 
gravitino mass of this order is favored from several things, i.e. the warm dark matter, the 
interesting region of the gaugino masses at the LHC, and the solution to the /i-problem. 
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A More on the extended SUSY breaking sector 

A.l Mass spectrum 
Scalar spectrum 

In order to analyze the spectrum of the scalar particles, we decompose the scalar compo- 



nents as 
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= -^{x<t> + iy^) ■ (43) 
In the foUowings, we suppress the index a of ^ and ys,m- 



^^lu this paper, we assume that all the parameters in the SUSY breaking sector are real valued, so 
that the CP-symmetry is not broken. However, notice that the CP violation in the hidden sector is 
insignificant for the CP violation of the SSM, because the phases of the gaugino masses are the same as 
(S*)* F , which can always be rotated away. 
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The squared mass matrices of x's are given by li 

^ / A^A^ + 2g\ {Sf -\'A {S) /V2 + 2^g\k {S) \ 

' 't4> [-X^A{S)/V2 + 2V2gj,A{S) 2X^A' + ^X^Sf + Agj,A^ J' ^ > 

for (xt, x^) and 

M^^)-=( -\\'A{S)/^ \ 

1^ -AA'A(5)/v^ (A'2 - A2)A2 + {Sf ) ' ^ ^ 

for [xsiXm)- On the other hand, the squared mass matrices of y's are given by, 



A^A^ -A^A {S) /v^ \ 



for (t/t, y<i,) and 



_^(,).^/ A-A^ -AA'A(5)/v^ \ 

^™ 1^ -AA'A(^)/v^ (A'2 + A2)A2 + iA2 (S)' j ' ^ ^ 

for ijjs.ym)- 

The eigen-modes of each mass matrices are given by. 



mi 



= ^ (^tr7W2^^(tr7W2)2_4detA^2^ . (48) 

Notice that the hghter mode of (j/t, Ucj,) is massless, which corresponds to the would-be 
Nambu-Goldstone mode of the gauged U{1) symmetry]^ 

Fermion mass spectrum 

The mass matrix of the fermion components of Sa and Ma is given by 

M^^^=i ° ^'^ A (49) 

The mass matrix of the fermion components of T, (p and the U{1) gaugino is given by, 

/ -AA ~iV2gx{S) 

-AA X{S)/V2 -i2gxA | • (50) 
V -tV2gx (S) -i2gxA 



The fermion component of S corresponds to the Goldstino. 



^® Note that the low-energy effective Kahler potential of is given by if ~ A^[((/) + (f)*)/\/2K 
'^gxVx)]^ + ■ ■ 'i where Vx is the gauge supermultiplet of U{l)x- 

^^There is no mixing between the R-axion and the would-be Nambu-Goldstone boson. 



23 



Vector boson mass 

The mass of the vector boson is given by, 

M'y = 2g\{{Sf + 2K^) . (51) 

The lightest particle 

As we have discussed, the model possesses a discrete Z4 symmetry which is effectively a 
Z2 symmetry under which Xm,s and ym,s are odd. Prom the above mass matrices, we find 
that the lightest, and hence stable, Z2 odd particle resides in Xm^g- In the foUowings, we 
name the lightest particle xl- 

In addition to the hghtest Z2 odd particle, the model possesses two massless scalars 
which correspond to the pseudo-fiat direction, a, and the R-axion, a, respectively. Besides, 
the model also has a massless fermion, the fermion component of S, which corresponds to 
the Goldstino. As we wiU discuss, the mass of the pseudo-flat direction is generated by the 
radiative corrections. The R-axion mass is also generated by the effects of the explicit R- 
symmetry breaking terms, e.g. the constant term in the superpotential generates the axion 
mass in supergravity. The Goldstino becomes massive by the super-Higgs mechanism of 
supergravity. 

Relevant interactions 

The relevant interaction terms to analyze the rehc density of the lightest Z2 charged 
particle is summarized below. 

First, the R-axion interactions only appear in the kinetic terms. In the basis we have 
defined above, the R-axion interactions come from the kinetic terms of S, T and Sa, 

The other important interaction term is 

^ = -^(^m + iym)Glpm + ^{Xm - iym)& , (53) 

where G denotes the Goldstino and ipm the fermion component of M^. 
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A. 2 Annihilation of the Hghtest Z4 charged {Z2 odd) field 

The annihilation modes of xl are xlXl GG, xlXl — >■ cc, and xlXl — >■ aa- In terms of 
xl, the relevant interaction terms can be rewritten by, 



{S) ml , A=^A^A(2A^+ (5)') 1 AU^-A^ {Sf 

2V2ml-mf'' 4 m% - ml "^^"'^'^ 16 {m% - ml){m% - X'^A^f 



2^2 



+ T77;^(^«)^^i + {cos9y{xLdyL - VlOxl) - sm9y{xLdyH - VhOxl)) 

4(d) v^W/ 

+^(aa)2 + ^/^^j + ^sin^^a;LG'(sin^/VL + cos^/V'/^) + ^-c. ■ (54) 

Here, O^^y are mixing angles of x and |/ components of Sa and M^, and 9f is the angle of 
the fermion components. The mixing angles are given by. 



-XX'A(S) ■ ^'-^'^ -AA-A(g> ■ ^^'''^^ -AA-A(g> ■'^^> 
A. 2.1 otl^^l ^ era 

The process, x^xl — > aa, proceeds via the t and li-channel exchanges of xl and xh as 
well as via the contact interaction. The amplitude of the t and li-channel xl exchange is 
given by. 



•M.2xL-*2a — W) 2 2 2 1 2 

2 \ mfr — mi \mi —t mi — u 



''H '"-LJ \<"-L '"-L 

. (56) 

2 {m% - ml) 

where we have neglected the mass of the fiaton and used t ~ m ~ —ml in the non- 
relativistic limit. The amplitude of the t and channel xh exchanges is given by, 

_ A^A- /2A^+(^\V 11^ 



16 \ mjj — ml ) \mjj — t mjj — 

l2 /o\2a2 I \2/c\2\2 / \/2A2 \ 



8 y mjj — ml J \ml-\-mjj^ 
The amplitude of the contact term interaction is given by. 



(57) 



M2.,-.2a - -T Ui-m? Ui-A'^A^ ■ ^^^^ 



''H~"''L/ \"''H 

Notice that the sum of all the above matrix elements vanishes for A = A'. 



25 



A. 2. 2 xlXl — )■ aa 



The process XlXl aa has four contributions, from the contact term interaction, the t 
and M-channel exchange of yi and the s-channel flaton exchange. The matrix element of 
the contact term interaction is given by, 

M.,.,^aa = - cos2 ^.^^ ^ '2 cos^ , (59) 

where we have neglected the mass of the R-axions in the final state. The matrix element 
of the t-channel exchange of y's is given by, 

_ cos el cos^ Oy / P3 ■ {-2pi + Pa) ■ {-2p2 + Pi) Pi ■ (-2pi + Pi) ps ■ {-2p2 + Ps) 



my — u 



4 cos 9^ cos^ 9y mj 



where we again neglected the mass of the R-axions. The contributions from the i/h 
exchanges are given by, 

4 cos^ 6'^, sin^ 6*^ 



Finally, the flaton exchange contribution is given by, 

mif — mf s 2 ml — mj 



A.2.3 xlXl GG 

The annihilation cross section into the gravitino is given by t and -u-channel exchange of 
Ma and Sa fermionso The matrix element is given by, 

2 /I 2/1 /-/ xit + ^\^'^ _/ xiu + ^^P / xl 

■^x^x^^GG = -r^^^ ^^^^^ ^yP^)^m V[pi)+U[p^)^jy^ V[pi) , (63) 

^ y nvi ' — t rrvl ' — u J 

where qt = Ps — Pi and Qu = P3 — P2- The heavier field exchange is then given by, 

,(/) ^ ^^(/) 



= X sin^ 9^ cos' 9, Up,)^-^±pi-v{p,) + u{p,)^-^^±p^v{p,) ] , (64) 



(60) 
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Here, we neglect the s-channel flaton decay which utilizes the higher dimensional operators. 
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Thus, the unpolarized squared amphtude is given by, 



\Mx,.,^ac\' - 2ASin^^.sin^g, / + 2ASin^ cos^ % 



-4A^ sin^ sin^ cos^ Of ^ '^^ . (65) 



Here, we have used the fact that the product of the two fermion masses is the negative 
valued since the determinant of the mass matrix in Eq. (H9l) is negative. 

A. 2. 4 Total cross section 

By adding up all the above modes, we obtain the total S'-wave cross section, 

|2 



1 1 \M[ 

A 

where a factor 1/2 represents the statistical factor of the final state particles. The thermal 
average can be trivially taken and the resultant cross section appearing in the Boltzmann 
equation is given by, 

1 1 \M.? 

A. 3 Decay of the charged particle 

As we have discussed in the paper, the Z2 symmetry is expected to be broken by the 
reduced Planck suppressed operators in Eq. ( I30ll . This operator generates the interaction 
terms of the Z2 odd particle such as, 

-> sin6'^A (5) - cos6'^A'A^ -^XLiHuHa)* - sin 6*^ ^^ XL^Hui^H^ + h.c. 

(68) 

Thus, the decay rate of the lightest Z2 odd particle is given by. 



Txr — — -smfe'a.A ^cost^2,'A , ^„ m^., H sm , mj., . (69) 

^ 47r \2 ^2 m^.^y M| ^ Svr Mj, ^ ^ 
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A. 4 The messenger interaction 

For the extended vector-like SUSY breaking model, we need at least two sets of the 
messengers, which couple to S±, 

kS m \ ( \ , ~_ \ / kS m \ ( '^2+ 



where we have rewritten S± in terms of the pseudo-flat direction S, and we have taken 
Thus, the Dirac-type fermion mass matrix is given by, 

\ m I 

for each \E'i's and \&2's. The squared mass matrix of the complex scalars is given by. 



/ k^{Sf + m^ km{S) kF \ 

km{S) m2 

kF k^Sf + m^ km{S) 

y km (S) j 



(72) 



The messenger sector has the U{l)di x U{\)a x U{l)d2 x f/(l)^2 global symmetries, 
which act on the down type and lepton type \E'i's and \l'2's, respectively. We have to break 
these symmetries explicitly to avoid the messenger dark matter overdosing the universe. 
For example, by introducing a small mixing such as 5m^i^2) we can break the symmetry 
down to U{l)d X U{l)i. Furthermore, we can completely break the symmetry by intro- 
ducing interactions such as e\I'il0ssM5ssM) where IOssm and Sgs^ are the SSM matter 
fields written in the 5'f/(5)GUT representations and e is a very small Yukawa coupling. 
What type of interactions are allowed depends on the precise R-charge assignment to the 
SSM matter fields, which we do not explicitly specify in this paper. 
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